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We consider the tunneling properties of a single fermionic impurity immersed in a Bose-Einstein 
condensate in a double-well potential. For strong boson-fermion interaction, we show the existence 
of a tunnel resonance where a large number of bosons and the fermion tunnel simultaneously. We 
give analytical expressions for the lineshape of the resonance using degenerate Brillouin-Wigner 
theory. We finally compute the time-dependent dynamics of the mixture. Using the fermionic tunnel 
resonances as beam splitter for wave-functions, we construct a Mach-Zehnder interferometer that 
allows complete population transfer from one well to the other by tilting the double-well potential 
and only taking into account the fermion's tunnel properties. 
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I. INTRODUCTION 

Bosc-Einstein condensates (BEC) have become a valu- 
able resource in today's research on many-body physics 
[l[ . Loaded into an optical lattice the low-energy regime 
realizes a Bose-Hubbard Hamiltonian [2[ whose param- 
eters can be tuned over a wide range by adjusting the 
optical lattice or engineering particle interactions via Fes- 
hbach resonances. 

A simpler, however, not less interesting variant of BEC 
in an optical lattice is obtained, when in the Mott- 
Insulator regime, the lattice is modulated by an addi- 
tional laser beam to create local double-well potentials 
at each lattice site. When tunneling between the local 
double-well potentials is negligible compared to tunnel- 
ing inside the double-well potential, the system is well- 
described by a two-site Hubbard Hamiltonian. BEC 
in double-well potentials have received considerable at- 
tention in past years Q. In Q, the dynamics are dis- 
cussed both on a mean-field level, where the non-linear 
Schrodinger equation becomes the discrete self-trapping 
equation Q, and in a quantum mechanically exact way. 
For large enough particle numbers, the two-site Hub- 
bard Hamiltonian can be approximated by the Josephson 
Hamiltonian [(| , which is the reason for these systems to 
be called Bose Josephson Junction (BJJ). A BEC in a 
double-well potential defines a representation of the ro- 
tation group and hence a pseudo-spin which can be uti- 
lized for quantum information processing and studies of 
decoherence Q . Other applications exploit the regime of 
weak tunneling, where the system behaves similarly to a 
quantum nanostructure in the sequential tunneling limit 

B|. 

A related strand of research is constituted by the in- 
vestigation of Bose-Fermi mixtures. These have been 
studied mostly on the mean-field level by solving non- 
linear Schrodinger equations for one or two wells (Tl| 



'Electronic address: Maximilian. Rinck@unibas.ch 



or by composite-fermion methods on an optical lattice 
[13 - [l4l |. Phase-diagrams have been computed for Bose- 
Fermi mixtures in three [l5| or one dimensions fl6l - [Tc| . 

In this paper, we want to study the dynamics of 
a Bose-Einstein condensate in a double-well potential 
when a single fermionic impurity is added to the sys- 
tem. The main question that comes into mind refers to 
the tunneling properties of the fermion when the double- 
well is tilted: Does the BEC leave the fermion's tunnel 
properties unaffected, or does the BEC expel the fermion 
to the other well against the potential gradient? Obvi- 
ously, the answer depends on the relative interactions 
between the two species. In our work, we consider the 
ground-state properties in the weak-tunneling limit. We 
discuss the different regimes in the parameter space de- 
fined by the particle interactions and the implications of 
large repulsion between the two species on the adiabatic 
and quasi-adiabatic dynamics. 

In Scction[Hj we introduce the Hamiltonian and its ba- 
sic properties. Then, in Section IlIIl we compute the ex- 
pectation value of the relative number operators indicat- 
ing that the ground state shows different phases defined 
by the repulsive forces. In particular, we find an avoided 
crossing between two states which are not connected di- 
rectly by the tunneling Hamiltonian. We calculate the 
splitting by an application of Brillouin-Wigner pertur- 
bation theory in Section IIVI In the last section of the 
paper, we consider time-dependent dynamics and show 
that by using the avoided crossings due to the tunnel- 
ing of the fermion as "beam splitter" , we can construct 
a Mach-Zehnder interferometer that allows us to trans- 
fer all population from one well into the other on a time 
scale that is only defined by the tunneling properties of 
the fermion. 



II. MODEL 

The starting point for our discussion is the standard 
Bose Josephson Junction [4| : a Bose-Einstein condensate 
is loaded into an optical double-well potential such that 
only the ground states of each well are occupied. Dcfin- 
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ing the relative number operator «b := j(tir — ^l) as 
the difference between the number of bosons in the right 
and the left well, the Hamiltonian is, up to a constant de- 
pending on the total number of particles N = jir + til, @ 

H B = -2en B + U B n 2 B - A B (b{b R + h.c.) . (1) 

The double- well potential can be tilted generating an en- 
ergy difference 2e between the two wells' ground states. 
The inter-particle repulsion of the bosons is of strength 
U B , and tunneling between the two wells occurs with am- 
plitude Ab- The operators hi, b\, annihilate and create 
a boson in the respective well. The two most promi- 
nent regimes of the Bose Joscphson Junction are the 
supcrfluid-likc regime U B <C Ab, where the particles are 
delocalized over the two wells, and the Mott-like regime 
Ab <C Ub, where the particles are localized. In the lat- 
ter case, tunneling is considered a perturbation to the 
Hamiltonian, which is diagonal in the relative number- 
state representation. Tilting the double-well potential by 
adjusting e, the bosons will tunnel into the other well one 
by one leadingto a staircase profile for the expectation 
value (ns) (e) H, 0- 

We now consider an additional single fermion, or, since 
in that case the particle statistics do not matter at all, an 
atom of a different species than the constituents of the 
condensate. Its dynamics are governed by the very same 
Hamiltonian without the repulsive term 

H F = -2en F - A F (dc n + h.c") , (2) 

with appropriately labeled constants and Cj, c\ being the 
fermion's annihilation and creation operators. We as- 
sume that the mutual interaction of both species is pro- 
portional to X] a =L r n B n F 0' wmcn in the relative num- 
ber representation reads 

Hb-f = 2U B - F n B n F (3) 

plus a constant. The full Hamiltonian of our system is 
thus 

H = Hb+H f +Hb-f. (4) 

The double-well potential in the two-mode approxima- 
tion defines a representation of the rotation group SU(2), 
hence a pseudo-spin on the Bloch-sphere whose length 
is the number of particles. The z-direction of this spin 
encodes the position information of the particles and is 
given by the relative number operator tjb/f with eigen- 
values m B /F [3- 

III. PHASE DIAGRAM 

As it is shown in [1, , by adjusting e adiabatically 
and thereby tilting the double-well potential, the Bose 
Joscphson Junction shows single-particle tunneling and 



a staircase profile of the expectation value of the rela- 
tive number operator. In our case, where an additional 
species, although only a single particle of it, is present 
in the system, we expect the same behavior in the tun- 
neling regime Ub , Ub-f 3> Ab /f '■ tilting the potential 
makes the particles tunnel from one well to the other. 
Since there are two species of atoms present in the po- 
tential, the obvious question to ask is, which species will 
tunnel? 

In general, transitions will be shown to be interaction- 
mcdiatcd, that is the relative magnitude of the repulsive 
interactions will determine the tunneling species. If there 
are no interactions as for instance in the case N = 1 , or at 
the exact threshold where the tunneling species changes, 
Ub = Ub-f, the profile of the transitions in (tib/f) is 
heavily influenced by the kinetic energy. 

In Fig. [1] we show numerical results for (n B / F ) as a 
function of inter-species interaction Ub-f and tilt e for 
two pairs of values of A j . The most obvious characteris- 
tic is the change of the tunneling species at the first res- 
onance for Ub-f = Ub- Close to this point, the fermionic 
expectation values show traces of attempted tunneling 
(a). For larger inter-species interaction, the bosonic ex- 
pectation values are shifted by one due to the presence 
of the fermion in the well with higher energy (b,d). For 
larger fermionic tunnel amplitude, the bosons show neg- 
ative compressibility kb = d(ns)/de close to zero tilt 
(d). In the following, we shall discuss and explain these 
phenomena in detail. 

Since we arc in the tunneling regime with well sep- 
arated resonances, we begin by restricting the bosonic 
Hilbcrt space at each resonance to the subspace spanned 
by the eigenstates of ub \m) = m \m), \m), \m + 1). The 
fermionic Hilbert space is also two dimensional, such that 
we have 

tt = (UBm 2 -2em -X m \ 

-A m U B (m + l) 2 -2e{m + l)J ' 

(5a) 

* = (-Af _ - A /) ■ ( 5b ) 

Hb-f = 2t/ B -F«B7iF (5c) 

in the respective basis. In the same basis, the relative 
number operators take the form 

" B = (o m + l)' " F= (7 (6) 

The off-diagonal elements of Hb are the tunnel ampli- 
tudes times the matrix element of the bosonic operators 

b a , bi, X m := AbsJ(±Nb + m + 1)(±/V B - m) [|. In the 
absence of tunneling, the Hamiltonian is already diagonal 
with energies 

E lm .±i } = U B m 2 ± Ub-ftu - 2e (m ± ^ . (7) 

If we place the system at a point, where tunneling of a 
particle should occur, we have either the states |m; ±5) 





and \m + 1;±|) degenerate or |m;±i) for a bosonic or 
a fermionic resonance, respectively. Equating the un- 
perturbed energies of the Hamiltonian, we find that the 
fermion can tunnel for e = itiUb-f, and the boson can 
tunnel for e = J7b(ti + ^) ± ^Ub-f depending on the 
position of the fermion. In particular, if we consider 
two values of the boson-fcrmion interaction, Ub-f,± with 
^{Ub-f,+ + Ub-f.-) = Ub, at U B -f = U B -f,+ , the 
state |m;— i) has the same energy as |m — 1; i) has at 
Ub-f = Ub-f,— This property is seen in Figs. QJb) and 
(d), where the bosonic resonances at equal distances left 
and right from Ub-f = Ub differ by one boson exactly. 

Approaching the first resonance of the system from 
large negative e, the ground state is \m = — N/2;— 5); 
if the fermion tunnels first, we have to equate its en- 
ergy with -E-at/2 ; i/2 otherwise with £7-jv72+i;-i/2- From 
the above formula for the energy, we see that the dif- 
ference between the positions of these resonances is 
£ F - £b = 5 [(1 - N)(U B -f - U B )}. Hence if U B -f > U B , 
the fermion will tunnel first and a boson otherwise. Re- 
peating this calculation for arbitrary m shows that this 
condition is independent of the resonance in question: 
for large enough inter-species repulsion, the fermion will 
tunnel first, otherwise last. 

This result is intuitively clear, as the condition states 
that for Ub-f > Ub, keeping the fermion together with 
the bosons costs more energy than keeping the additional 
boson. Hence the fermion will be expelled to the other 
well. Noting that such behavior occurs for finite e, the 
fermion will move to the potential well with higher en- 
ergy, hence against the potential gradient. 

Directly at the degenerate point Ub-f = Ub, it can- 
not be decided from the atomic interactions alone, which 
particle will tunnel first. In this situation, we have to 
include the tunnel Hamiltonian, as now the kinetic en- 
ergy of the particles will decide. At e = ttiUb = itlUb-f, 
the unperturbed Hamiltonian's ground-state is four-fold 
degenerate: 

\™> - !; 2) ' \ m >~\) > K \) A m + !; -\) 

all have the same energy. Including the tunneling the 
Hamiltonian in this basis reads 
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FIG. 1: (Color online) Number expectation for a single 
fermion and iVbosons = 11- Left: (jif), right: (ub)- The 
parameters are Af = 0.01 and Ab = 0.05 in the upper plots 
and exchanged for the lower plots. The color code is Black- 
Red- Yellow for Negative-Zero-Positive. 



For the first resonance from positive/negative e, the state 
|m ± 1; =pi) drops out, and the expression simplifies to 



(«F> = ± 
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(10a) 
(10b) 



Note that these are discrete values evaluated at e = 
ttiUb-f and are not continuous in m. 

From these expressions, we directly infer that at the 
first resonance, a fast tunneling species, i.e., one with 
large A^, will hamper the other and restrict its num- 
ber expectation value to the asymptotic value of ±i or 
±A/2, respectively. 



plus a constant. This Hamiltonian has a bi-quadratic 
characteristic polynomial, which can be solved explicitly, 
leading to the ground-state expectation values 



(n F ) 



Aim 



^/Ap + 4A 4 B m 2 + AgAp [-4m 2 + N(N + 2)] 

(9a) 

(ub) = m - (jif). (9b) 



By considering the limit of large and small tunnel- 
ing amplitudes, we can significantly simplify the expres- 
sions (|9l^l9b"|) . Also, as for Ub-f < Ub only single-particle 
processes are present, the situation at Ub-f = Ub can 
be used as a good approximation for the ground-state 
properties in that regime. For Ub-f > Ub we observe 
processes where a larger number of particles tunnels si- 
multaneously. This will be discussed in Section IIV1 In 
the present case, the asymptotic expectation values at 
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for Ab/Af — > 
- — S for A B = A F 

Vl+AT(iV+2) 

isgn(m) for A B /A F — > 



oo 



and 



(«b) ->• < m I 1 - 



for Ab/Af 
for A R = A r 







m — o sgn(m) for Ab/Af — > oo. 

For Ab = Ap, shown in Fig. [2] (top panels), the fermionic 
expectation value shows steps linear in the resonance 
number m. For Ab < Ap, Fig. [5] (middle panels), the 
fermionic expectation shows oscillations about (rip) = 0. 
In the other case, Ab > Ap, shown in Fig. [2] (lower 
panels), we see that the presence of the fermion and in 
particular its tunnel amplitude do influence the bosons 
in such a way that for Ab ^S> Ap, i.e., a strongly local- 
ized fermion, the relative number expectation value ap- 
proaches m — g sgn(m) which is the same as in the case 
Ub-f < Ub- Also the expectation- value of the fermionic 
relative number operator approximates a step function. 

In the case of low fermionic tunnel amplitude, the res- 
onance at zero tilt e = 0, with m = will also show neg- 
ative compressibility «b = d(ns)/d£ < 0, because in the 
limit A F = 0, \im m ^ (nB) = \ and lim^o («b) 
see Fig. [3] A quantitative result is given in 
< e < Ub-f < U B - 



i 
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for 



IV. ZERO-BIAS SPIN FLIP 

Due to the symmetries of the Hamiltonian without tun- 
neling, Ab/ f = 0, at e = 0, the states \niB',rn-p) and 
|— TOb;— nip) are degenerate. In particular, this holds 
for the ground state of the non-interacting system. In 
general such a degeneracy for the ground state always 
occurs at a resonance. In our system, however, the de- 
generate states do not only differ by a single particle 
having changed its position, but rather correspond to 
the exchange of positions for two species and, depending 
on the inter-species repulsion Ub-f, a larger number of 
bosons. Since, as we have mentioned above, the relative 
number operators are equivalent to the z-components of 
the pseudo-spin defined by the double- well potential, this 
transition corresponds to a flip of the combined pseudo- 
spin of bosons and fermion. 



A. Higher-order degenerate perturbation theory 

In contrast to the single-particle resonance, whose 
physics is that of an avoided crossing and degenerate 
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FIG. 2: (Color online) Occupation number difference (m) 
(bosons left/fermions right) for a mixture of 10 bosons and 
one fermion. The plots are for different values of the Bose- 
Fermi repulsion Ub-f G {0.75, 1, 1.25} in red/blue/black. The 
tunneling amplitudes are Ab = Af = 0.05 (top), Ab = 
0.01, A F = 0.05 (middle), and A B = 0.05, A F = 0.01 (bot- 
tom). 
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FIG. 3: (Color online) Negative compressibility kb = 
d(ns)/de due to low fermionic tunnel amplitude Af = 10~ 4 , 
A B = 0.05, and [/ B -f /Ub = 0.9. 



perturbation theory, the zero-bias spin flip is not easily 
amenable to degenerate perturbation theory, because the 
involved states are not directly connected by the tunnel- 
ing term in the Hamiltonian. Although a numerical treat- 
ment of the problem is straightforwardly implemented, 
perturbation theory allows for a simple analytical ap- 
proach to the resonance and provides a very good ap- 
proximation for the lineshape. 

For general considerations, let \m) and |— m) be eigen- 
states of the unperturbed Hamiltonian with energies 
E± m — £o i e such that they are degenerate for e = 0. 
Let £y he a, perturbation, with a scalar £ to keep track 
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of orders of magnitude, and (— m | V | m) — 0. Assume 
the shortest chain of matrix elements of V that connects 
the two degenerate states via intermediate states \rij) has 
length k. Brillouin-Wigner perturbation theory to lowest 
order in £ yields for the splitting (l9j 



1 



E„ 



E, 



"»ni,n2 



1 



-^712 ^Ti 



(n) 



The lineshape of the expectation value of the relative 
occupation number n = m(\m) (m\ — \—m) (— m|) re- 
stricted to the degenerate subspace as a function of e 

is thus to lowest order in £, (n)(e) ss |m|e/ -i/ \ A£" 2 n + e 2 . 



In the regime with C/b-f > Ub, the spin-flip transition is 
|— to; ^) i-> |to; — i). To construct the chain of intermedi- 
ate states connecting the two degenerate states by single- 
particle processes, we climb up the angular-momentum 
ladder in to. At a certain to = too, the fermion jump 
is included. Then, due to the different repulsion ener- 
gies, the energy-denominators are different and the chain 
naturally splits into two products. Consider the energy 
differences 



G< X (n) :=E n 



E 



= {n 2 - ni 2 )U B + (n + m)U B -F (12a) 
= {n 2 - m 2 )U B - (n - m)U B -F (12b) 



B. Application to the Bose— Fermi mixture 

Let us apply this reasoning to the spin flip at zero 
bias of a Bosc-Fcrmi mixture in a double- well potential. 



for bosonic transitions left (G^ 1 ) and right (G>*) of mo. 
The fermionic jump is given by G> 1 (too), which leads to 
the result 



J 



A£ = 2A F A fe £ 



\k= — m / m =—rn 



G<(n) ] [ G>(r 



(13) 



The fermionic jump is included as the first factor of the 
last product. Here, we use the convention that a product 
without factors, eg., Yim w ith n < to, is unity. 

As an example, we consider the smallest non-trivial 
mixture: N = 2. For Ub-f > U B , the spin-flip transition 
at zero bias is | — 1; 1/2) M> [1; — 1/2}, which are connected 
by three paths 




The level-splitting is thus 

AE = 4A R A F 

Ub-Ub-f 



U B - 



^ II;-: 



1 

U B - Ub-f 



(14) 



(15) 



and, since to = 1, as a function of the tilt e, (jib)(£) = 
e/^Je 2 + |A£; 2 . In Fig. H we show the numerical data 

for kb = d(riB)/de as well as the line shape computed 
with Brillouin-Wigner perturbation theory. 



Of course, the result is not absolutely accurate, as in 
the tails of the resonance, we do not recover the asymp- 
totic states with well-defined occupation numbers, but 



rather the exact eigenstates of the full Hamiltonian. If 
we look at the numbers of Fig.[U Brillouin-Wigner theory 
gives AE = 3 • 1CP 3 , whereas the numerically evaluated 
splitting is approximately 2.88 • 10~ 3 . 



V. LANDAU-ZENER DYNAMICS 



In the previous sections, we have focused on the adia- 
batic dynamics of population transfer from one well into 
the other by increasing the tilt e of the potential. In a 
realistic scenario, however, we will always adjust the tilt 
within finite time, hence with finite velocity de/dt. This 
means that we have to take into account Landau-Zencr 
physics of quasi-adiabatic transitions (20l - [24| . Depend- 
ing on the velocity and the splitting of the states at a 
resonance, the population transfer is thus heavily influ- 
enced... 

In the regime Ub-f > Ub of our model, the zero-bias 
spin-flip resonance is so narrow that it could be very dif- 
ficult to traverse it adiabatically as the necessary veloci- 
ties, depending on the splitting, are very small. In [25| . 
such an effect is also seen and successfully exploited to 
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FIG. 4: (Color online) The Brillouin-Wigner result (dots) 
and the numerically obtained lineshape (solid line) for the 
compressibility kb = d(na)/d£ of the zero-bias spin flip. Note 
that the energy scales have to be separated by about two 
orders of magnitude J/U ~ 10 -2 . Here Ub = 1, Ub-f = 2 , 
A B = 0.1, and A F = 0.05, such that AE = 6A|A F . 



achieve the population transfer sought for. Let us choose 
a velocity a for the change of e(t) = at in time that al- 
lows to pass the single-particle resonances adiabatically. 
Starting at infinite negative time, the fcrmion tunnels 
from left to right at finite negative e as it is expected for 
Ub-f > Ub- The spin-flip resonance is, however, passed 
completely diabatically, such that all particles stay where 
they are and do not exchange places. At this point the 
system is no longer in its ground state but in the first 
excited state. Further tilting of the potential causes the 
bosons to tunnel to the lower-lying well step by step. Fi- 
nally, the system adiabatically passes an avoided crossing 
of two excited states, where again, the fermion changes 
sides and tunnels back into the left, i.e., higher- lying well, 
as is seen in Fig. [SJa). There it will stay until other pro- 
cesses than those describe by our Hamiltonian cause it 
to tunnel again. The diabatic traversing of the spin-flip 
resonance thus causes the transition of the system from 
its ground state at t — >• — oo to the first excited state at 
t — >• oo. Lowering the velocity allows us to pass to the 
adiabatic regime also for the central resonance, but the 
more bosons we have in the condensate, the slower we 
would have to adjust e(t). 

In order to facilitate adiabatic population transfer from 
one well to the other, we make use of a constructive inter- 
ference effect. The spectrum of the Bose-Fermi mixture, 
as it is shown in Fig. [5] is symmetric about e = 0. If we 
adjust the tunnel amplitudes such that Ap < Ab, we can 
control the Landau-Zener physics of the fermionic jump 
and still pass all other bosonic resonances almost adia- 
batically. This allows us to restrict the problem to only 
consider the ground state \g) and the first excited state 
|e) of the Hamiltonian as depicted in Fig. [5] for the case 
of TV = 3. Passing the first avoided crossing by adjusting 
the tilt e{t) amounts to a unitary transformation of the 



asymptotic states 

jj _ ( cos © sm @\ 
I— sinO cos0y 

with the angle depending on the velocity a and 
the splitting AE of the two states cos 2 = 1 — 

cxp (~^f ^cT~) 01 • ^ ne spin-flip transition is very nar- 
row, such that we traverse it diabatically, hence exchang- 
ing ground and excited state, which corresponds to the 
unitary transformation a x . The second avoided crossing 
is identical to the first; however, it is passed in opposite 
direction, whence the invoked unitary is U'. Since both 
arms of the interferometer are identical — they can be 
mapped by e n- — e onto each other, the accumulated dy- 
namical phase is the same and amounts to a global phase 
factor, which is of no importance for our purpose. 

With these assumptions, the constructed Mach- 
Zchnder interferometer is the map 

Ua x U f = a x cos(20) - a z sin(20) 

For = 7r/4, U<t x W = — er z : the ground state of the 
full Hamiltonian will again be mapped onto the ground 
state and population can be transfered completely. Con- 
sidering the avoided crossings as a beam-splitter for an 
incoming ground-state wave function 23\, this choice of 
amounts to a velocity a for which the Landau-Zener 
transition splits the wave-function exactly in two halves 

\g)^^(\g) + \e)). 

In the numerical data in Fig.[5{b), the achieved popu- 
lation transfer is almost perfect. The expectation values 
of the relative number operators show some oscillations 
in the fermionic part, which are due to the coherent su- 
perposition created by the Landau-Zener transition. The 
remainder of the profile is, however, well in accordance 
with the adiabatic picture for atom counting, only that, 
since the wave function is split by the avoided crossing, 
we do not count full but half atoms, as the change in (m) 
is 1/2 instead of the expected 1. 

The important advantage of this approach over a com- 
pletely adiabatic transfer is that with this interferome- 
ter, we only need to adjust the Landau-Zener transition 
of the single fermionic resonance. This, however, is in- 
dependent of the number of bosons in the system, and 
we can use a much higher velocity than in the purely 
adiabatic regime where the Landau-Zener physics of the 
spin- flip transition are taken fully into account. 

VI. CONCLUSION 

In this paper, we have investigated a Bose-Fermi mix- 
ture in a double-well potential with the restriction that 
while the number of bosons is arbitrary, the number of 
fermions is fixed at one. We thus discussed the influence 
of a single fermionic impurity on the ground-state prop- 
erties of the Bose Josephson Junction, when the potential 
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FIG. 5: (Color online) Expectation values (m) of the relative number operator for fermions (green/solid) and bosons 
(blue/dashed) in a numerical Landau-Zener experiment on the Bose-Fermi mixture. Ab = 0.1, Af = 0.05, Ub = 2, Ub-f = 3.5, 
(a) a = 0.001, (b) a « 0.011. 
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e 

FIG. 6: (Color online) Spectrum of the Bose-Fermi mixture 
and the path that is traversed by a linear scan of e. The 
branches that define the two equivalent arms of the Mach- 
Zender interference are shown in red/dashed and green/dash- 
dotted, black/solid denotes the in- and outgoing beam. 

was allowed to be tilted. It has turned out that rather 
than the bosons, it is the fermion that is affected most by 
the boson-fermion interactions. We have separated the 
dynamics into two regimes: one, where the bosons and 
the fermion live side-by-side, and one, where the fermion 
is expelled from the BEC towards the higher-lying well. 
In this regime, we have found a zero-bias spin-flip, where 
the fermion and a large number of bosons change places. 
The physics of this process, in particular the level split- 
ting and the lineshape of the observed resonance in the 
relative particle number, are well accessible by Brillouin- 
Wigner perturbation theory. Since many particles are in- 
volved in this transition, complete population transfer of 
all particles from the left to the right well upon tilting 



the potential is hardly possible anymore, when we assume 
the realistic case of slow but not infinitely slow adjust- 
ment of the potential tilt e. Instead, we have shown, how 
the fermionic resonance in the regime [/b-f > Ub can be 
employed as a beam-splitter of a Mach-Zchndcr inter- 
ferometer to achieve population transfer at much higher 
tilting speeds than would be necessary to traverse the 
zero-bias spin-flip resonance adiabatically. 

There are several direction in which we want to extend 
our investigation. In the case of very large TV, the Gross- 
Pitaevskii equation provides a much simpler description 
of the BEC than can be achieved by numerical diagonal- 
ization of the two-site Bose-Hubbard Hamiltonian. We 
will therefore ask for the lowest-order corrections of the 
dynamics of a Bose Josephson Junction with large N in 
the presence of a single fermionic impurity. The trade- 
off for the reduction to the Gross-Pitaevskii equation is, 
however, the extension to the study of nonlinear eigen- 
value problems. 

The second direction refers to the number of fermions 
in the double-well potential. Experiments of two 
fermions with opposite spin in such a system without 
a BEC have been conducted by Trotzky et al. [26| . The 
orbital wave function of two fermions with opposite spin 
in a double-well potential is decomposed into a singlet 
and three triplet states with respect to the pseudo-spin 
defined by the double-well potential. Interestingly, the 
singlet does not couple to the bosons at all. On the con- 
trary, the state |f, -J-), that is a configuration with definite 
position of the spin-up fermion in the left well and the 
spin-down particle in the right well, is a superposition of 
the singlet and the J z = triplet, which, however, docs 
couple to the BEC. In this configuration we expect the 
time-dependent dynamics to show very interesting phe- 
nomena, which could be observed experimentally. 

This work has been supported by the Swiss SNF and 
the NCCR Nanoscience. 
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